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ABSTRACT 

Accretion disks that become gravitationally unstable can fragment into stellar or sub-stellar com- 
panions. The formation and survival of these fragments depends on the precarious balance between 
self-gravity, internal pressure, tidal shearing, and rotation. Disk fragmentation depends on two key 
factors (1) whether the disk can get to the fragmentation boundary of Q=l, and (2) whether fragments 
can survive for many orbital periods. Previous work suggests that to reach Q=l, and have fragments 
survive, a disk must cool on an orbital timescale. Here we show that disks heated primarily by ex- 
ternal irradiation always satisfy the standard cooling time criterion. Thus even though irradiation 
heats disks, and makes them more stable in general, once they reach the fragmentation boundary, 
they fragment more easily. We derive a new cooling criterion that determines fragment survival, and 
calculate a pressure modified Hill radius, which sets the maximum size of pressure-supported objects 
in a Keplerian disk. We conclude that fragmentation in protostellar disks might occur at slightly 
smaller radii than previously thought, and recommend tests for future simulations that will better 
predict the outcome of fragmentation in real disks. 



1. INTRODUCTION 

The past two decades have witnessed growing inter- 
est in the conditions under which gravitationally unsta- 
ble disks fragment to produce bound objects. These 
fragments may grow into massive stars or black holes 
within AGN disks, or small stellar and substellar com- 
panions in protostellar disks. In either context, we can 
divide the process into two steps. First, the disk must 
become Hnearly gravitationally unstable, and secondly, 
once formed, the fragments must be able to survive the 
various disruptive forces in the disk. 

The first step, linear gravitational instability (GI), 
requires that the disk be Toomre unstable: Q = 
CsVL/ttGY. < Qo ^ 1 (jToomrd [1961 ISafronovl [l960h . 

where Cg = y^kT/ja is the isothermal sound speed of the 
gas with mean particle weight ja and temperature T, G is 
the gravitational constant, S is the disk surface density, 
k is the Boltzmann constant, and ft is the orbital angular 
frequency. 

The second stage is more uncertain. Three primary 
forces compete to destroy newly born fragments: internal 
gas pressure, tidal forces from the central object (shear- 
ing), and internal rotation. Fragments can also destroy 
each other through mutual interactions. 

Whether the disk can be driven to the fragmentation 
boundary, which we call Q = Qo^ and whether fragments 
can survive, must be considered separately. Both re- 
quirements place limits on the cooling properties of disks. 
Moreover, the answer to both of these questions depends 
on whether the disk is actively heated (i.e. by dissipa- 
tion of accretion energy through an effective turbulent a 
viscosity), or passively heated (i.e. be external irradia- 
tion). Both processes operate in most disks, but it mat- 
ters which is the dominant energy source. Here we con- 
sider the passive (irradiated) case in detail, as it has been 
neglected in previous theoretical studies. Throughout 
the text we use the term viscosity to describe any process 



that contributes to both transport of angular momentum 
and local energy dissipation. 

The role of cooling in fragmentation has been discussed 
in the liter ature in two way s, which we dub the Viscous 
Criterion (Pringle" 1981; Gammie '2001) and the Colli- 
sional Criterion (Shlosman & Begelman 1989). The Vis- 
cous Criterion describes whether an actively-heated disk 
can rid itself of the energy generated by Gl-driven ac- 
cretion on the disk's dynamical timescale, while the Col- 
lisional Criterion considers whether newborn fragments 
will collide with each other on this timescale. While the 
Viscous Criterion addresses the first question: (1) Can a 
disk be driven to Q = Qo so that initial collapse begins, 
the Collisional Criterion addresses the second: (2) Can 
a collapsing fragment survive? 

Here we provide a more complete picture of the nec- 
essary disk conditions for fragmentation and fragment 
survival. We begin with question (1), emphasizing the 
crucial role of disk history in determining whether a disk 
can be driven to collapse. The Viscous Criterion answers 
this question solely for viscously- heated disks, and we 
provide a complementary criterion for irradiated disks, 
which we show automatically satisfy the Viscous Crite- 
rion. We then turn to question (2). We find that the 
Collisional Criterion does not fully capture the role of 
cooling in fragment survival, because it does not account 
for the contraction of fragments. We propose a more 
general cooling criterion which governs fragment survival 
in all disks. This "Stalling Criterion" first determines 
whether fragments collapse on the free-fall timescale, or 
if they quickly become pressure supported. If a frag- 
ment becomes pressure supported, stalling collapse, it 
must be smaller than a pressure-modified Hill radius to 
avoid tidal disruption. Fragments that contract to small 
enough sizes before stalling can survive while slowly un- 
dergoing Kelvin-Helmholz contraction. 

Fragment survival depends both on the cooling time 
and the gas equation of state. The Stalling Criterion is 



applicable in both viscously- heated and irradiated disks. 
Figure [T] illustrates how the different cooling criteria we 
describe determine whether disks fragment, and whether 
their fragments survive. 

We begin with brief definitions of the disk cooling time 
and related terms in Section O We then distinguish 
quantitatively between disks heated via external irradi- 
ation or internal viscous dissipation in Section 13.11 and 
describe the run-up to instability in each case (Sections 
13.21 and 13. 3p . We introduce our new cooling criterion 
in Section |4l We show that for arbitrary cooling time, 
the Stalling Criterion sets a critical size scale (i?H+p) be- 
low which a fragment can survive if pressure supported. 
Whether or not a fragment can reach this critical radius 
depends on the disk equation of state. We review other 
means for disruption of fragments in Section [5l 

We discuss the implications of these processes for pro- 
tostellar and AGN disks in Section [6l then comment 
briefiy on recent simulations of viscously-heated disks 
(Section [7j). Finally, we provide suggestions for future 
simulations of irradiated disks and conclude in Section 

E 

2. THE DEFINITION OF THE COOLING TIME 

For clarity we begin with a few definitions. In general, 
we consider a disk with cooling time 
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where Tc is defined by the following term in an energy 
equation: 

du u 

(2) 

cool '^c 



dt 



/3 is a constant, and u is the internal energy per unit 
mass. The use of the /3 cooling formalism is motivated 
by the fact that the orbital time is the important dynam- 
ical timescale, and allows for a cancellation in the energy 
equation. The result will be a dimensionless cooling cri- 
terion. 

In reality, Tc is a function of the physical properties 
of a fflven disk, and thus varies in time and location. 
iKratter et al.l (|2QlQb[ ) calculate the cooling time for a 
thermal perturbation from the midplane to b43* 
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where tr is the optical depth due to the Rosseland mean 
opacity, tr = hzIl/2 , 7 is the adiabatic index, a is the 
Stefan-Boltzmann constant, and we have assumed that 
the disk is optically thick, which is valid for Q = 1 disks. 
We use Equation (0) to translate between dimensionless 
cooling times and physical disk conditions, however we 
quote our analytic expressions in terms of Tc for general- 
ity and convenience. 

With these definitions in hand, we now describe when 
astrophysically relevant disks are susceptible to fragmen- 
tation. 

^ In strongly viscously-heated disks, an order unity correction 
factor must be applied to the coefficient of Equation ^ if the 
opacity of disk material depends on temperature (see K ratter et al.l 
[2010b ). Including this factor does not change our results. 



3. DRIVING DISKS TO LINEAR INSTABILITY 

The required initial conditions for creating gravitation- 
ally unstable disks depend on the dominant energy source 
in the disk. As alluded to above, passively heated disks 
behave differently from actively heated ones. First we 
clarify this distinction, and then review the required con- 
ditions for disks in each class to become linearly unstable 
toGI. 

3.1. Irradiated vs Viscously Heated Disks 

What does it mean for a disk to be "dominated" by 
irradiation? Qualitatively it means that the primary 
source of the disk's thermal energy is external irradia- 
tion rather than dissipation of accretion energy (which 
we compute as an artificial viscosity). In other words, 
the temperature to which the disk midplane is heated by 
external irradiation is higher than the equilibrium tem- 
perature that the disk could achieve by balancing accre- 
tional heating and cooling via thermal emission. As a 
result, irradiated disks will be roughly vertically isother- 
mal. We shall show an example of this later in Figures [2] 
and [31 

We choose this terminology over the more common, 
but less informative, distinction between passive (irradi- 
ated), and active (viscously heated), primarily because 
irradiated disks can have ongoing accretion - they are 
not entirely "passive." 

Quantitatively, the boundary between an irradiated 
and viscously heated disk can be defined by balancing 
radiative cooling from the midplane with heating via dis- 
sipation of accretion energy. We quantify this dissipation 
in the standard manner, by employing an effective vis- 
cosity, v = ac^/Q, described by the parameter a < 1. 
The midplane temperature, T, generated by dissipation 
of energy due to accretion of a mass fiux 



M = SttSz. 



is given by: 
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in the optically thick case, relevant for Q = 1 disks. In- 
serting Equation (J5j) into Equation (J3j) yields 
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Equation ([6|) expresses the conditions for thermal equi- 
librium of a disk heated only by viscous dissipation in 
terms of the cooling time as defined in Equation dU . Ir- 
radiated disks have cooling times less than this critical 
value: 

^cool < — -^~ , (7) 

iza 7 — 1 

because they are heated to higher temperatures by ir- 
radiation than by dissipation of accretion energy. This 
critical disk temperature, in terms of density and a is: 
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If stellar irradiation maintains the disk (or some region 
of it) above the temperature in Equation (J8]) , it is irra- 
diation dominated. Typically the outer regions of disks 



Q=l? 



Fragment Survives? 




criterion satisfied? (^ 



Do spiral arms iceep Q>1? 



Are fragments disrupted? 



Fig. 1. — Flow chart illustrating the relevant physics and criteria for determining disk fragmentation. First, the disk much reach Q = 1. 
In viscously heated disks this is determine d by the Viscous Criterion f ^3.2|) . and in irradiated disks by whether global modes can process 
mass as fast as it in falls on to the disk f ^3.3|) . If Q — ^ 1, then one must determine if fragments survive. Collapsing in free fall likely 
guarantees survival f^4.1D . Fragments which do not collapse quickly must avoid collisional f ^5.1|) or tidal disruption by shrinking below a 
critical size, i?H+p (" ^4.21) . 



are irradiation dominated, while in their inner regions 
accretional heating is more important than irradiation. 
This radial dependence is due to the liberation of more 
energy as material in the disk is lowered deeper into the 
potential well of the central object. 

To evaluate whether a disk in steady-state satisfies 
Equation ([8]) at a particular disk location, one must set 
a to the value required to process the incoming material 
either from larger radii in the disk or the background 
environment (e.g. a protostellar core). Of course, if the 
a required to process the material is larger than that 
provided by any means of angular momentum transport, 
a steady state is not possible. We note that over time, 
a particular disk radius may transition from irradiation- 
dominated to viscously heated or vice versa as the disk 
evolves. 

The boundary between viscously heated and irradiated 
disks is independent of the mechanism for angular mo- 
mentum transport, so long as it can be quantified by an 
a. In this work we are concerned with transport due to 
GI. Unlike trans port due to the magnet o-rotational in- 
stability (MRI) (;Balbus h Hawlevlll99 ^. which is often 
modeled using a single value for a, Gl-induced transport 
cannot be described by a fixed a. Instead, a is a func- 



tion of disk parameters (in particular Q) and should be 
thought of as a response to t he self-gravitating state of 
the disk (jLin fc Pringlel[l987[ ). We return to this point 
in Section [3^21 

We now answer the first question posed above, and 
in Figure [T] how do disks reach the boundary of linear 
instability? 

3.2. Driving a Viscously-Heated Disk to Instability - the 
Viscous Criterion 

In order to fragment, a disk must first reach Q = Qo ^ 
1, allowing linear collapse to proceed. Though dramatic 
events may sometimes drive disks into gravitationally un- 
stable configurations on short timescales, typical disks 
likely evolve on long enough timescales that Q = 1 must 
be approached gradually if it is approached at all. 

First, consider disks which are heated only through ac- 
cretion energy. In the absence of any external heat source 
(e.g. irradiation ) , a disk can only be dr iven to fragment 
by rapid infall (jMatzner fc Levinll2QQ5f ). In the case of 
a protostellar disk, for example, this infall comes from 
the protostellar cloud and its rate is set by conditions 
external to the disko We discuss in the next section how 

^ Evolution to Q = Qq might also occur in protostellar disks 



irradiation changes this picture, and in the appendix we 
derive exphcitly the evolution of Q with dM/dt. 

As long as infall rates are not too large, a Q ^ 1 disk 
heated solely by viscous dissipation relaxes into an equi- 
librium state known as gravitoturbulence. Gravitotur- 
bulence results from a feedback loop which may be un- 
derstood as follows. When Q ~ 1, if radiative cooling 
exceeds viscous heating, the disk will cool and get closer 
to Q = 1. As it gets closer to Q = 1, the onset of 
gravitational collapse leads to turbulence and accretion, 
which heats the disk and generates extra pressure sup- 
port to counteract gravitational collapse. If the heating 
rate were to exceed the cooling rate such that Q ^ 1, 
then the a generated by GI would decrease, the heating 
rate would decline, and the disk would again cool down. 
This feedback loop is gravitoturbulence. 

In gravitoturbulent equilibrium, viscous heating and 
radiative cooling balance such that 



the Viscous Criterion: 
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(|Gammiell2QQll ). Here, /3, defined in Equation ([T]) , repre- 
sents the cooling rate in terms of orbital periods. One 
may derive Equation (|9]) by equating the cooling rate per 
unit mass u/tc (Equation[2]) with the viscous heating rate 
per unit mass, (9/4)ac^r^ (c.f. Equation [5l remembering 
to include two disk faces). We have used -u = 0^/(7 — 1), 
appropriate for an ideal gasQ. 

However, the equilibrium reflected in Equation (J9]) can- 
not be achieved for arbitrarily high mass-infall rates. For 
Gl-driven accretion, a is a function of disk properties (in 
particular of Q), dictated by the strength of self-gravity 
in the disk. Numerical simulations demonstrate that 
there exists a maximum transport rate provided by grav- 
itational instability (Gammie 2001) and hence a maxi- 
mum a < 1, which we call its saturation value agat- 

The precise value of Ofgat as a function of disk prop- 
erti es, and the function a(Q) remain un certain (though 
see, iKratter et all 120081 : IVorobvovl I2010D . We empha- 
size that a is not a free parameter, and neither is the 
temperature of the disk in equilibrium. Whether or not 
a disk can restablize itself at ever higher M depends on 
both a{Q) and on the temperature dependence of the 
opacity. Unless hz ex T^^, increasing M always drives 
Q down. Depending on the functional form of ^(Q), a 
very steep temperature dependence may also cause Q to 
decline with increasing M. The disks that w e consider 
here typically have k, oc T-^ (jSemenov et al.l 12003). but 
we explore this dependence fully in Appendix A. Thus in 
general, we expect rapid accretion to drive disks toward 
instability, and in some cases fragmentation. 

The distinction between a disk that can maintain grav- 
itoturbulence and one that fragments is encompassed in 



after infall from the protostellar cloud ceases if the efficiency of 
surface accretio n due to t he MRI varies as a function of distance 
from the star ( Zhu et all '2010. Perez-Becker & Chiang, in prep). 
The following still applies, however, with surface accretion from 
outer regions of the disk adopting the roll of "infall." 

^ The extra factor of 7 in the denominator of Equation (|10|) . 
included for consistency with Gammie (2001), comes from that 
paper's use of the adiabatic sound speed rather than the isothermal 
sound speed in the definition of Q. 



When Equation (p!Q|) is satisfied, cooling exceeds viscous 
heating and the disk is unable to reverse this imbalance 
by increasing its accretion rate, and hence its heating 
rate. Thus material is being added into the disk faster 
than it can be drained to smaller radii, which causes S 
to increase proportionally faster than T^/^, Q to decline, 
and fragmentation to proceed. We note that accretional 
energy must be deposited in the disk as material falls 
toward the star, independent of the physical processes 
involved in dissipation of that energy. Hence, the Viscous 
Criterion is robust unless dissipation of accretion energy 
is dramatically non-local. 

Note that most simulations of gravitoturbulence in- 
volve isolated disks, i.e. not undergoing infall. These 
disks fragment because they are set up out of equilibrium 
(jLodat o & Rice 2004, 2005). Cooling exceeds heating by 
fiat, and thus the disk moves towards Q = 1. In reality, 
the cooling rate is set by the thermal physics of the disk, 
and is not an independent parameter. 

Of course, as we addressed above, many disks have 
two heating sources, accretion energy and external irra- 
diation. In fact, the Viscous Criterion is nearly identi- 
cal to the boundary between these two regimes in Equa- 
tion (O, modulo different numerical coefficients. The dif- 
ference in coefficients derives from the fact that here we 
are concerned with the timescale for the disk to radiate 
away its internal energy, and therefore cool off. To dis- 
tinguish between viscous and irradiated disks we instead 
consider whether a temperature perturbation in the mid- 
pl ane heats the disk or i s radiated away (see the appendix 
of IKratter et al.ll2010b[ ). A disk satisfying the Viscous 
Criterion can avoid fragmentation only if it is irradiation 
dominated, because the external heat source alters the 
feedback loop - the temperature, and thus cooling rate 
- are no longer controlled by dissipation, but instead by 
the radiation field. We address how disks heated by ir- 
radiation are driven towards instability below. 

3.3. Driving an Irradiated Disk to Instability 

What is required for an irradiated disk to be driven to 
Q = Qo? Again, infalling material is required to drive 
the disk unstable. However, we suggest that for very 
low mass disks, the critical infall rate may be lower than 
previously thought. 

Unlike their viscously heated counterparts, irradiated 
disks will not become gravitoturbulent as infall adds 
mass and decreases Q. Gravitoturbulence relies on a 
feedback loop between Q and disk temperature, but the 
temperature of an irradiation-dominated disk is indepen- 
dent of Q, modulo small perturbations. For simplic- 
ity, for the remainder of this section, we only consider 
disks which remain irradiation dominated at all accre- 
tion rates. 

To understand why irradiated disks cannot maintain 
gravitoturbulence, consider the following hypothetical 
disk, with Q = 1 around a 1.5Mq star. Imagine that 
as the disk approaches Q = 1, the effective a generated 
by Gl-driven accretion increases from 0.1 to Ofgat = 1.0. 
For an order of magnitude increase in a, the equilibrium 



temperature for a disk heated through viscous dissipa- 
tion rises from 5 — 20K at 100 AU. However, the equi- 
hbrium temperature of the disk including external irra- 
diation from a young A star increases from 47 — 50K. 
Thus the fractional contribution of GI to the total pres- 
sure support of the disk (ex T) is negligible. Now, we 
add more and more material to the disk. Even if the 
disk processes material at this maximum rate, and dis- 
sipates energy at this maximum rate, the disk cannot 
reheat itself via turbulence. External irradiation, which 
initially prevents the disk from collapsing, cannot adjust 
the disk temperature to avoid the Q = 1 threshold. Once 
Q crosses this threshold, fragmentation will proceed. 

However, irradiation-dominated disks are not entirely 
unaffected by self-gravity as infalling material drives 
them toward Q = 1. Irradiated disks (or indeed 
disks with an isothermal equation of state, which be- 
have similarly) respond to decreasin g Q by generating 
gravitational torque s in spiral arms (jBolev et al.l l2007c 
iKratter et al.l 1201 Oaf ). These torques can, in turn, drive 
rapid accretion through the disk. 

The properties of spiral arm transpo rt dep end on the 
ratio of the disk mass (Adams et al. Il989f ). Mdisk = 
TrrQ^^S, to the central mass, M*. Here rout is the outer 
radius of the disk, where most of the mass resides. For a 
Q = 1 disk, Mdisk /^* = H/rout^ where H = Cs/^ is the 
disk scale height, so the disk-to-central object mass ratio 
is equivalent to the disk aspect ratio. Spiral arms that 
form in thin and thick disks have different properties. 

Numerical simulations indicate that high-mass (thick) 



disks can process niass infall rates up to M 
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(IKratter et al.ll2010at[Qffner et al.l l2010V while at higher 
rates they fragment. Lower infall rates do not generate 
fragmentation because these disks become globally grav- 
itationally unstable; the resulting spiral arms transport 
material efficiently, keeping Q > Qq. In the language of 
Section [321 ^ ^ ^l/^ corresponds to driving the disk 
with an infall rate having a ^ 1. Thus thick, irradiated 
disks, like viscously- heated disks, must be forced with 
large infall rates in order to fragment. 

However, if the disk is low in mass (thin), has inefficient 
transport by spiral arms, and is very near Q = 1, a mod- 
est infall rate could tip it over the fragmentation thresh- 
old faster than it can process even the modest amount of 
mass it receives. For this to occur, other mechanisms of 
angular momentum transport such as the MRI must not 
be able to process the infalling material. Whether low- 
mass (thin), irradiated disks can avoid effective angular 
momentum transport by spiral arms remains uncertain. 
Simulations suggest that strong gl obal spiral arm trans- 
port may require Mdisk /M* > 0.1 (JLodato fc Ricel [2004). 

Since the actual criterion for global mode transport is 
unknown, it is difficult to provide a precise criterion: fu- 
ture numerical work can delineate the exact conditions 
for this regime, and indeed whether it exists at all. This 
work will need to be conducted with great care; N-body 
and SPH simulations which show spontaneous instabil- 
ity at Q > 1 may suffer from Swing amplification of 
Poisson noise (D'Onghia et al 2011, in prep) . Current 
state of the art simulations (|Bolev et al.ll2007l : iCai et all 
12008^ do not yet explore the full parameter space re- 
quired to answer these que stions. For earlier attempts t o 
quantify this behavior see iLaughlin fc Rozvczkal () 19961 ) ; 



iLaughlin fc Korchadnl (jl996f ): iLaughlin et al.l (jl997aD 

Can nature produce a Q ~ 1 disk that remains 
irradiation-dominated, unaffected by the MRI, and low 
enough in mass to avoid global instability? We show in 
^that it is difficult, but not impossible to achieve Q = 1 
and maintain a sufficiently low mass disk to suppress spi- 
ral arm formation. 

Thus we have answered the first of our two questions 
and shown how both classes of disks can reach the insta- 
bility threshold. In summary, to fragment, viscously- 
heated disks must be driven by mass infall at a rate 
corresponding to a > o^sat- Disks that are irradiation- 
dominated (and thus satisfy the Viscous Criterion) might 
be driven to fragment at lower accretion rates, particu- 
larly when such disks have low masses relative to their 
host stars (i.e., are thin). This possibility is interest- 
ing because such disks could generate fragments that 
remain near their initial masses rather than experienc- 
ing substantial growth. Numerical simulations of global 
mode formation in low- mass, irradiated disks will be re- 
quired to evaluate this possibility. Other angular mo- 
mentum transport mechanisms such as the MRI will pre- 
vent fragmentation from occurring for arbitrarily small 
infall rates. 

Can fragments survive once the instability threshold is 
reached? We now consider the cooling criteria related to 
fragment survival: the Stalling Criterion, and the Colli- 
sional Criterion. 

4. DETERMINING FRAGMENT SURVIVAL 

Long term fragment survival requires that the frag- 
ment's own self-gravity wins out over internal pressure, 
rotation, tidal forces from the central object, and interac- 
tions with fellow fragments. We consider the importance 
of the first three forces through the Stalling Criterion, 
then argue that in all regimes of astrophysical interest, 
fragments which satisfy the Stalling Criterion are also 
immune to collisional disruption. 

The Stalling Criterion concerns the evolution of an 
already- formed fragment — in other words, it constrains 
the non-linear outcome of successful linear collapse. If 
a disk can reach the fragmentation threshold, two pos- 
sible fates befall a newly formed fragment. If the frag- 
ment cools on its own free-fall timescale, it can contract 
quickly to a length-scale that is small compared to its 
Hill radius. Conversely, if cooling is slow compared to 
the free-fall timescale, then the fragment's collapse will 
be stalled by pressure support, and it will begin to col- 
lapse on a modified Kelvin-Helmholtz timescale, of order 
its cooling time. 

We expect a fragment that collapses on the free-fall 
timescale to survive, although we discuss other mecha- 
nisms for disruption in ^ Will a fragment whose col- 
lapse is stalled by pressure support on the disk dynamical 
time survive? It depends on the radius at which pressure 
halts collapse. To survive, a pressure supported frag- 
ment must be smaller than a critical size scale, which 
we call the pressure- modified Hill radius, Ru-\-p- Once 
a fragment has shrunk to a size smaller than i^H+p, it 
can survive regardless of its cooling rate (modulo other 
mechanisms for disruption discussed in ^. 

Shrinking quickly to a fraction of a Hill radius may 
prove impossible for some slowly cooling fragments, de- 
pending on the gas equation of state and on the initial 



density perturbation seeding collapse. We calc ulate the 
requirements for collapse on a free-fall time in §4.H de- 
rive i?H+p in E21 and estimate the initial conditions 
required for fragment survival in §4.31 

We comment that the Stalling Criterion has been over- 
looked primarily because gravitoturbulence is thought to 
generate a large a^at- If Q^sat is large, a disk which cools 
quickly enough to be either irradiation dominated or sat- 
isfy the Viscous Criterion (see Equations 171 and [TQ|) will 
usually cool quickly enough to ensure free-fall collapse of 
the fragment. However the physics controlling fragment 
survival is distinct from that controlling the approach to 
the (5 = 1 boundary, and for smaller Ofgat, or disks with 
sharp opacity jumps, both criteria should be considered. 

4.1. Survival via Free-Fail Collapse 

Consider a disk that has arrived near enough to Q = 1 
to produce perturbations that can seed collapse. In an 
isothermal disk, once a perturbation exists with Q < 1, 
collapse will occur on the free fall time scale. For non- 
isothermal collapse we must examine the energy budget 
of the collapsing region as self-gravity tries to make it 
contract. 

We neglect internal rotation in the following analysis 
as it is difficult to estimate its evolution in time analyti- 
cally. We comment on the importance of rotation in ^ 
By neglecting rotation we provide optimistic limits for 
collapse and thus fragment survival. 

When Q ^ 1 and linear collapse has begun, self-gravity, 
pressure, and tidal forces from the central object exert 
comparable forces on the fragment, but self-gravity is 
currently winning. As collapse proceeds and the radius, 
X of the fragment gets smaller, the ratio of the self- 
gravitational force to the tidal force increases as x~^, 
assuming that the enclosed mass in the fragment stays 
constant. However, the ratio of self-gravitational force to 
pressure force may increase or decrease. We define the 
virial ratio 5 in a clump with density p: 
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where Menc = 4/37rpx^, and P is the internal pressure. 
We approximate the pressure gradient as P/x since the 
fragment is only one pressure scale height in extent 
If B increases with time, then self-gravity will always 
dominate and free fall collapse is ensured. If B decreases, 
pressure will eventually stall collapse. 

For small enough /3 (fast cooling), B initially increases 
as collapse proceeds. In other words. 
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where the negative sign comes from the fragment radius, 
X, decreasing as the object collapses. Fragments collapse 
in free fall when (3 is less than a critical value, which we 
now calculate. 

Holding fragment mass, Menc, fixed and substituting 
in for an ideal gas with pressure P = (7 — l)pii. Equation 

^ Including tidal gravity in B causes a small change in the derived 
coefficient for the free-fall regime. We exclude it here for simplicity. 



12]) yields the following inequality: 

__d_ f l_ 
dx \xu 



>0. 



(13) 



To proceed we need an expression for the spatial 
derivative of u: 



du 
dx 



P_dp 



u 

Tr. 
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dt 
dx 



where the time derivative of p is given by 



dp 
'dt 



p dx 
X dt 



(14) 



(15) 



The coefficient 3 in Equation ([15]) arises from the spatial 
derivative of the enclosed mass and equals the dimension 
of collapse. 

The three terms in our expression for du/dx (Equation 
[T4]) come from PdV work, cooling parametrized by the 
cooling time Tc (c.f. Equation [2]) , and viscous heating. 
Balancing viscous heating against cooling gives rise to the 
Viscous Criterion. In order to isolate the effect of PdV 
work after fragmentation has occurred, we shall here- 
after drop the last, viscous term, assuming that turbu- 
lent heating is inefficient within a fragment. This choice 
differentiates fragments from turbulent eddies. 

For convenience we define the collapse timescale of the 
fragment: 

/ dlnx \ ~ , , 

icon - - (^ j (16) 

Equation ([13]) may now be written in the following simple 
form: 

tcoii>rc(37-4) (17) 

Equation ([TT]) confirms that the object cannot collapse on 
a timescale faster than its own cooling time. We recover 
the standard result from stellar astrophysics that if 7 < 
4/3, thermal pressure cannot prevent collapse. Also, an 
object cannot collapse on a timescale faster than its own 
free-fall time, 



te 



v/2^ 



(18) 



where to derive the final expression, we have assumed 
that p = T,/2H with H = Cg/O. equal to the disk 
scale height. Hence, the fragment collapses on whichever 
timescale is longer: Tc (37 — 4) or tff . 

Comparing these timescales, we find that free fall col- 
lapse requires: 

As expected, pressure is negligible when the cooling time 
is of order an orbital time. For Q = 1 the limits on j3 for 
free fall collapse are: 



/?^13,7: 



5/3 
7/5. 



(20) 
(21) 



These critical values of j3 will be compared to the actual 
radiative cooling time of fragments in different regimes 
ingi 



This calculation involves order unity coefficients which 
remain uncertain. Our primary goal is to establish the 
relevant physical parameters that must be evaluated to 
determine if fragments survive; in reality each value de- 
pends on the non-linear evolution which follows linear 
instability. High resolution hydrodynamic simulations 
will provide more insight than an analytic consideration 
of a range of fragment structures, densities and cooling 
laws. Note that our use of j3 does not imply that the 
fragment cooling time must be equal to the background 
disk cooling time, although we expect the two to be quite 
similar just after the onset of fragmentation, which is the 
timescale addressed by this calculation. 

4.2. A Pressure Modified Hill Radius 

Now we consider the opposite regime, in which cooling 
is slow, and the fragment undergoes Kelvin-Helmholtz 
contraction rather than free-fall collapse. We calculate 
the critical size, i^H+p, at which a pressure supported 
fragment can withstand tidal shear. In the absence of 
pressure support, this radius would just be the Hill radius 
(or Roche lobe). The addition of pressure support is a 
small, though potentially significant, modification to this 
limit. Pressure partially negates the force of self-gravity, 
leaving the object more susceptible to disruption by tidal 
forces. 

In binary stellar systems, modification of the Roche 
lobe due to pressure support is small because the scale 
height at the surface of the star is much smaller than 
the size of the star itself. For fragments, this is not the 
case — the pressure scale length of the gas is comparable 
to the size of the fragment and pressure should be con- 
sidered. Thus we are justified in our approximation of 
the pressure gradient here as P/x. We calculate i^H+p 
for a compressible gas along the line of sight to the star, 
the direction in which tidal gravity is most disruptive. 

Consider a fragment that has reached the Kelvin- 
Helmholtz contraction regime. The system is evolving on 
a cooling time, which we assume is substantially longer 
than a dynamical time. This implies that self-gravity, 
pressure, and tidal forces are roughly in balance: 
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(22) 



Because the fragment is contracting slowly, it remains in 
approximate equilibrium: 



d fGMencP 
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(23) 



and x{t) is an undetermined function describing contrac- 
tion. For simplicity, we rewrite this expression in terms 
of the dimensionless variable B: 



d_ 
dx 



B 



\Rh) 



= 



(24) 



where Rh = a (Menc/3M*)^/^ is the Hih Radius. 

We plug Equations ([14]), ([15]), and ([22]) into Equation 
([23]) to calculate the contraction time as a function of Tc 
and fragment size. Equation (|24|) reduces to: 



(-fefl(¥--^)«fe)'-<-' 



Substituting in for different values of 7 we find: 
4^-1 



^coll = < 



where 
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,7 = 5/3 
,7 = 7/5 



(26) 



(27) 



We can find the critical value of y for which the clump 
will keep contracting by finding where tcoii, given in 
Equation (|26|) . changes sign: a negative contraction time 
signals expansion. This provides an estimate of the max- 
imum size fragment that can be pressure supported but 
not torn apart by tides: 
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1/3 



RHn = 7/5 



(28) 
(29) 



4.3. The Stalling Radius vs Ru-\-p 

To interpret the severity of the hurdle embodied in 
Equations (|28|) and (|29|) . we must relate i?H+p to the 
initial size and mass of the fragment. A debate regarding 
the appropriate choi ce of these number s is ongoing in the 
literature (see, e.g.. iBolev et al.l [20101 ). The choices we 
make here are not the only sensible possibilities. We 
take the diameter of the collapsing region, I^coib to be 
half of the Toomre wavelength. At, at Q = Qo ^ 1, 
envisioning that only half of the wavelength participates 
in forming an overdensity. Since At = 27rH, ^coii = 
7tH. We expect this length scale to be at a minimum 2H 
for collapse to be roughly symmetric in the vertical and 
planar directions^ The radius of the collapsing region 
^coii = ^coii/2 = At/4 = (7r/2)JY, and the mass enclosed 
in the cylinder circumscribed by this radius is: 

Mfrag = Sq^i^(At/4)^ (30) 

The Hill radius for a fragment with this mass is: 
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(31) 



(32) 



Thus At/4 cannot be an estimate of the initial radial 
extent of the fragment for any extended length of time. 
Rather, At/4 likely represents the initi al feeding zone for 
the fragment, not the fragment itself (jRafikovll2005l ). A 
collapsing fragment will only appear as a bound object 
when X < Rh- 

This interpretation is not an artifact of our choice of 
rcoii- In a Q = 1 disk, the only length scale on which one 

^ In reality, the scale height and thus wavelength may depend 
on the struc ture of the spiral arm in which the fragment forms 
(|Bolev et al.l 2010 ). 



can create an object of mass Sttx^ that lives within its 
own Hill radius is smaller then the disk scale height: 



Rh,sc = TTT^ 



H 
3Q 



(33) 



Creating such a fragment is inconsistent as it implies that 
the vertical and radial extent of the fragment differ by a 
factor of three. Apparently, one cannot view fragments 
as chunks of the disk of radius x extracted as with a 
cookie cutter from a sheet of dough. Perturbations must 
grow from the Q = 1 state, and contract. 

Can such a contracting fragment achieve the radius 
required by Equation (|28|) or ([29]) before pressure stalls 
collapse? In general, answering this question requires 
following the non-linear evolution of the fragment, in- 
corporating changes in the cooling time throughout this 
process. Full consideration of this evolution is beyond 
the scope of the current work. However, we obtain some 
guidance by considering the collapse of a fragment in the 
adiabatic (i.e. infinitely slow cooling) limit. 

Imagine that at the time of its initial linear collapse, a 
fragment has virial ratio B = Bq > 1 but —dB/dx < so 
that pressure is growing faster than self-gravity. Assum- 
ing that rotational support does not become important 
(see Section [5]), collapse will stall when B ^ 1. The frag- 
ment will survive if tidal forces are not important at this 
stalling radius (i.e. if the stalling radius is smaller than 
i?H+p)- For an adiabatic gas, P oc p^, so 



B oc x(3^-^) 



(34) 



where we have assumed a constant Menc so that p ex 
x~^. Hence, the fragment stalls when its radius has been 
reduced by a factor of approximately Bq^ for 7 = 5/3 or 
Bq^ for 7 = 7/5. 

For rcoii = 7rH/2 and Qo = 1, ^0 = 7r/2. The frag- 
ment will initially start collapsing. In a 7 = 5/3 gas, the 
fragment only collapses to ~60% of its original size in 
the adiabatic limit (cooling would allow additional col- 
lapse). However, for 7 = 7/5, the collapse factor is nearly 
10. In fact, to contract from At/4 ~ 1.7 Rh to OARh 
(c.f. Equation [29|) requires only a ~ 34% enhancement of 
gravity over pressure {Bq ~ 1.34) in the limit of infinite 
cooling time. This large factor suggests that even if the 
cooling time of the fragment becomes long, if its equa- 
tion of state is soft, it can survive. This further suggests 
that fragment survival is trivial in disks that have both 
soft equations of state (including 7 = 7/5, appropriate 
for molecular gas) and strong external irradiation. 

That the 7 = 5/3 case is less promising for fragmen- 
tation is neither surprising, nor as significant. We have 
only considered the adiabatic limit. Nevertheless we infer 
that fragmentation is much more likely to be inhibited by 
the Stalling Criterion for stiffer equations of state. Note 
that the slow cooling limit differs from the adiabatic case 
in that the former is not dissipationless. A truly adia- 
batic fragment could rebound after contraction. We use 
this limit merely to show that slowly cooling fragments 
can contract below the stalling radius for soft equations 
of state. Again, Equation ([3l|) reveals that fragmenta- 
tion is difficult to inhibit in radiation pressure dominated 
disks with 7 = 4/3. T his result is consistent with recent 
simulations by iJiang fc Goodman , (2011,) . 



5. OTHER DISRUPTION MECHANISMS FOR FRAGMENTS 

Thus far we have focused on thermal barriers to 
fragmentation. We now address dynamical and long 
timescale processes which may disrupt fragments. In 
particular we address the second classic cooling criterion 
mentioned in the introduction, the Collisional Criterion. 

5.1. The Collisional Criterion 

iShlosman fc BegelmanI ()1989f ) argue that slowly cool- 
ing fragments remain susceptible to disruption by neigh- 
boring fragments if the disk breaks up into objects sepa- 
rated by of order their own size, as expected from linear 
theory. Since the time to shear across a Hill radius at the 
Keplerian shearing velocity is of order an orbital period, 
fragments must contract on roughly the orbital timescale 
in order to avoid collisions. For a Q = 1 disk, the Ke- 
plerian shear energy at a separation of the Hill radius 
is of order the binding energy of a fragment, implying 
that encounters are likely d isruptive. More generally, 
IShlosman fc BegelmanI (|1989f ) suggest 



/? 



< 



(35) 



where Sq=i is the critical surface density for fragmen- 
tation. This relation derives from the proportionality 
between S and the wavelengths of growing modes in an 
unstable disk. Higher values of S provide a wider range 
of growing modes. 

This argument ignores an important effect, which is 
that for soft equations of state (e.g. 7 = 7/5), fragments 
with p > 1 can still collapse on an orbital timescale. 
Moreover, for fragments stalled at i^H+p, the binding 
energy should be larger than than the collisional energy 
by a factor of two, so fragments may avoid collisional 
disruption (and may not collide at all) for soft equations 
of state. 

Global simulations suggest that in systems with disk- 
to-central mass ratios ^ 10%, fragme nts form in rela- 
tive isolation within spiral arms (.Lodato fc Ricd I2QQ5I : 
iBolev et al.l I2QQ7I : iKratter et all |2QlQd), reducing: the 
overall frequ ency o f collis ions. Finally 'Goo dman fc Taiil 
(2004) and iLevinI ()2007l ) have argued that collisions 
might cause growth rather than disruption. These 
collisions have be en observed in some global models 
(jBolev et al.l[20Toh . 

5.2. Rotation 

The most likely candidate for disruption that we have 
neglected is rotation. Gauging the initial angular mo- 
mentum of a fragment and predicting its evolution is 
challenging. If angular momentum is conserved in a 
collapsing fragment, the object will eventually become 
rotationally supported. Dekel, et al., (in prep) report 
that fragments in isothermal disks collapse by a factor 
of 32/7r ^ 10 in radiu s befor e becoming rotationally sup- 
ported . iBolev etldl (|2010f ) found that in SPH simula- 
tions with radiative cooling, the rotational support was 
only 20% of the gravitational binding energy. 

In reality, subdisks eventually form around fragments, 
transporting away some fraction of the internal angular 
momentum. Since the disks are small, the viscous time 
is relatively short, about ^disk/^disk allowing them to 



drain faster than the parent disk. Recent numerical work 
also suggests that gravitational torques can likely prevent 
the object from rotating faster than 50% of break up 
(jLin et alJl2Qllf ). It is also possible that in the process 
of angular momentum transport, some mass is carried 
away. We make no attempt to account for this. 

If the above estimates are correct, then rotational 
stalling likely occurs at radii smaller than i?H+p, imply- 
ing that rotational support does not inhibit fragmenta- 
tion. However, more numerical work is needed to confirm 
this conclusion. 

5.3. Migration 

Even if a fragment can survive at its current loca- 
tion, migration may bring it closer to the host star, and 
shrink Ru+p below the size of the fragment. Because 
migration-induced disruption occurs on timescales longer 
than a few orbital times, even initially free-falling frag- 
ments might in principle suffer from late disruption due 
to this process, depending on the extent of the migration 
(jVorobyov & Basu 2005; Boley et al. 2010). 

5.4. Interactions with Spiral Arms 

Close interactions with spiral arms generate a tidal 
force roughly equivalent to that from the central star, so 
collapsing in free fall or shrinking to Ru+p should protect 
fragments from this means of disruption. In the vicinity 
of an arm of width Aarm, the ratio of tidal gravity from 
the arm to tidal gravity from the star is approximately 

(SarmAarm/^*)(^/'^arm) ^ (^arm /Aarm )/(7r(5 arm)- 
6. IMPLICATIONS FOR FRAGMENTATION 

We have argued that if a fragment can shed its angular 
momentum quickly enough to avoid rotational support, 
fragmentation may be generically possible in an irradi- 
ated disk with 7 = 7/5, appropriate for the outer regions 
of protostellar disks. This conclusion applies so long as Q 
can reach its critical value, Qo ^ 1. Here, we consider the 
conditions under which disks can be driven to Q = Qo, 
while remaining irradiation dominated. We comment on 
the implications for fragmentation in protostellar disks 
(Section 16. ip . AGN disks (Section 16. 2p are less affected 
by these considerations. 

6.1. Protostellar Disks 

We consider the case of a fragmenting disk surrounding 
an A-star. This choice is motivated by both theory and 
observations which suggest that A-stars may have com- 
panion s consistent wi th form ation by gravitational insta - 
bility (Kratter et al.' ^2008,, i20TQbl : iHinklev et aT][20Toh . 
We calculate the equilibrium disk temperature for a 
Q = Qo = 1 disk, including both viscous heating and 
irradiation from the star. 

We begin by setting a = Ofgat, its saturation value 
for the gravitoturbulent regime, and we identify the 
irradiation-dominated regions within the disk given this 
choice. Because agat is not well-determined, we con- 
sider a range of possible values. Because gravitoturbu- 
lence cannot support a > agat, disk regions which are 
irradiation-dominated for a = a^at will be irradiation- 
dominated in all Q ~ 1 disks. In contrast, regions of 
the disk which are viscously heated in this limit are not 
necessarily viscously heated at lower mass infall rates. 



Minf, onto the disk, so after considering this limit, we 
turn to the lower values of a which can be generated by 

low Minf. 

We use the stellar irradiation model of 
iChiang fc GoldreichI (J1997I ). which determines the 
incident flux on to the disk as a function of radiufl 
The disk midplane temperature is found by solving the 
following equation: 



crT^ = ^TrFj, + o-7]^rad 



(36) 



where the energy flux due to viscous dissipation in the 
midplane is 



Stt 



(37) 



the mass accretion rate through the disk when Q = 1 is 
M ^ 3acl/G, and 

'^\V'r-'/' . (38) 
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In the following plots, we use 
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where fi = 2.33mH and mn is the mas s of a hydro- 
gen atom. From the tables accompanving Tandin et aTl 
(2009), a pre-main-sequence A star has intrinsic luminos- 
ity L^ ^ 2.5L0 for ages of 1-3 Myr, with Lq equal to 
the current luminosity of the Sun. Infall is much more 
likely to drive disks unstable at earlier times, and at 3-8 
xlO^ years, the intrinsic luminosity of a 1.5 Mq star is 
approximately 7 Lq. Accretion luminosity will exceed 
the intrinsic stellar luminosity for accretion rates onto 
the star > lO~^M0/yr. The stellar accretion rate may 
or may not match the mass infall rate Minf- We choose 
L^ = BLq as an intermediate value. We use the cooling 
time defined in Equation (|3]) . 

To derive the most optimistic limits on fragmentation 
at small masses, we imagine that dust grains have grown 
to ~ 300/im so that the opaci ty k. ^ 0.24 cm^^~^ is 
constant with temperature (see lKratter et al.ll2010bl for 
a discussion of this choice and why it provides optimistic 
limits at small fragment masses). Given this disk model, 
we find the characteristic fragment mass STrr^^jj for the 
disk where fragments can free fall collapse. Again to 
generate the smallest pos sible fragment masses , we use 
rcoii = 7tH/2 (c.f. Section 101 iBolev et al.ll2010f ). 

Figure [2] shows the temperature profiles, critical radii 
where the cooling criteria are satisfied, and initial frag- 
ment masses for a disk with a^at = 0.05, We show the 
corresponding temperatures, radii, and fragment mass 
scales for Ofgat = 0.005 and 0.1 in Figure [H For Ofgat > 
0.05, which is expected for Gl-driven turbulence, the lo- 
cation at which the disk becomes irradiation dominated 
is effectively the same as the fragmentation boundary 
typically calculated from the Viscous Criterion. 

For lower values of agat fragmentation is pushed in- 
wards, and there are disk locations that remain irradia- 

^ We note that this is not entirely self-consistent for a disk 
heated in part by accretion, but the effect is subtle. 
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Fig. 2. — Top: Equilibrium and viscous disk temperature for 
an A star as a function of radius for a disk with gravitoturbu- 
lent a = asat = 0.05, for L* = 5Lq. The vertical labelled lines 
show the oft used critical /3 = 3 (green), the dividing line be- 
tween viscous, and irradiation dominated heating (blue), and the 
regime where free-fall collapse (FFC) is ensured for 7 = 7/5 (grey). 
Bottom: corresponding minimum initial fragment masses (black) 
where fragmentation is allowed and corresponding enclosed disk- 
star-mass ratio (grey) implied by the constant Q model. Note that 
fragments will continue to grow — low initial masses do not guar- 
antee the formation of planets over brown dwarfs. 







Radius (AU) 

Fig. 3. — Similar to Figure [2] with temperatures, critical radii, 
and fragment masses shown for disks with a = asat = 0.005 (solid) 
and 0.1 (dashed). Note there is only one line labelled FFC, corre- 
sponding to the a = 0.005 case because when a = 0.1, /3 < 13.5 
at all radii shown. In the large a case, the viscous regime extends 
to larger distances. The thin portion of the fragment mass curve 
for a = 0.005 indicates the points at which free fall collapse is not 
guaranteed, but the disk remains irradiation dominated. Fragmen- 
tation in these locations is more uncertain. Again, these represent 
initial minimum masses; fragments will grow significantly towards 
the gap-opening and isolation mass (,Kratter et al.M2010b. ') 



tion dominated, but where cooling prohibits free-fall col- 
lapse. Fragmentation is more uncertain in these regimes, 
but plausible if Ofgat ever takes such low values. 

Recall that the critical value agat depends on the 
transport mechanism, so in irradiation-dominated re- 
gions where global modes control transport, the satura- 
tion value for angular momentum transport will be dif- 
ferent from the value set by gr avitoturbulen ce and likely 
depends on disk mass (Laughhn et al."1997b). Low mass, 
irradiation dominated disks with low infall rates might 
be able reach Qo as discussed in Section [331 In this case, 
fragmentation could be pushed significantly closer to the 
star. 

To determine whether this low disk mass, low accretion 
rate regime might be relevant, we calculate the equilib- 
rium temperature for a Q = 1 disk given a constant 
accretion rate (not fixed agat), and find the irradiation 
dominated regions of the disk. In Figure S] we show the 
temperatures, disk masses, fragment masses, and corre- 
sponding values of a for three different accretion rates. 
We must compare these with aeff generated by the MRI, 
as noted below. 

As discussed in Section 13.31 Q = Qo is only likely 
to be achieved for irradiated disks with Md/M^ = 
H/r < 0.1, evaluated at r = rout, where this criti- 
cal value of Q.l remains substantially uncertain. Using 
the .Chiang fc Goldreid] |l997.) irradiation temperature 
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where we have assumed that the disk is optically thick 
at the outer edge when Q ^ 1. From the pre-main- 
sequence evolutionary tracks of iLandi n et al. (2009|), 
stars with M* < 2Mq have intrinsic luminosities L^ « 
3.5L0(M*/M0)i-^ at an age of 3-8 xlO^ years and 
L* ^ IALq{M^/MqY-^ at 1-3 Myr. These scalings 
suggest that at a given age, the value of H/r decreases 
somewhat with increasing stellar mass, implying that spi- 
ral arm-driven transport may be more effective for disks 
around G stars than for those around A stars, particu- 
larly when low accretion rates do not generate significant 
stellar accretion luminosity. 

Figure 2] illustrates that disks can plausibly be low 
enough in mass to avoid global instabilities, but require 
more angular momentum transport than provided by the 
MRI. The quantity Ofeff is averaged over the entire disk 
column — in dead zones experiencing only surface accre- 
tion, the aeff which can be supported by the MRI can be 
quite low ( ^ammia 1996.). In this case, the disk could be 
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Fig. 4. — Top: Black curves show the equihbrium disk tempera- 
ture for three different infall rates, Mini- The three blue, vertical 
lines indicate the radius outside of which the disks are irradiation 
dominated, while the grey vertical line (labeled FFC) shows where 
free fall collapse is ensured (it remains the same for all accretion 
rates). Middle: Disk and fragment masses as a function of radius. 
As in Figure O the thin lines indicate where the disk is irradiation 
dominated but free fall collapse is inhibited. Bottom: Effective a 
values corresponding to the chosen temperature and accretion rate. 
Even low values of a^f^ could exceed those f or layered MRI driven 
accretion above a dead zone (|Gammielll996l ). 

driven unstable while accreting at a lower rate than re- 
quired in the viscous regime. For such low rates to drive 
fragmentation, the disk must already be near Q = 1 such 
that low accretion rates need not significantly alter the 
disk mass. 

Since the exact criterion for global mode transport is 
unknown (see Section [X3|) , it is difficult to evaluate the 
importance of this regime. We suggest that future sim- 
ulations search for the maximum mass at which Q ^ 1, 
irradiated or isothermal disks do not show spiral arm- 
driven transport. 

We have shown regimes in which planetary mass frag- 
ments might form and survive. We emphasize that newly 
born fragments must still be prevented from accreting 
surro unding disk material if they are to remain low in 
mass (|Kratter et al.|[2QTQb[ ). The isolation mass for such 
fragments remains in the brown-dwarf to stellar mass 
regime. 



6.2. AGN and Galactic Disks 

Contrary to the protostellar case, the question of inter- 
est in galactic scale disks is "can fragmentation be pre- 
vented?" rather than "can fragmentation be relevant for 
planets?" Relaxing the cooling criterion in the strongly 
irradiated portions of disks is of little import as these 
regions already satisfy the more restrictive V iscous Cri- 
terion for large as at ^ even at a fraction of a pc (|Goodmanl 
l2QQ3l : iLevinI [20071 ). If galactic disks are irradiation dom- 
inated at radii as small 0.01 — 0.1 pc, fragments still col- 
lapse in free fall for modest black hole masses (lO^M©). 
The inner 0.01 — 0.1 pc of hotter, more massive AGN 
disks {Mbh = 10^ — 10^) are likely dominated by radi- 
ation pressure at these radii, shifting to 7 = 4/3, where 
free fall collapse in ensured . Rece nt numerical simula- 
tions by I Jiang fc GoodmanI (J2011l) are consistent with 
this analysis: they find fragmentation at long cooling 
times when radiation pressure dominates. 

7. RELEVANCE OF THE STALLING CRITERION IN 
VISCOUSLY HEATED DISKS 



In a recent suite of SPH simulations, 'Meru fc Batd 
()2010f ) find that the critical value of (3 given by the Vis- 
cous Criterion increases with increasing resolution, with 
no convergence in sight. The authors suggest a possible 
breakdown of the Viscous Criterion. A complete break- 
down of the Viscous Criterion would imply that dissipa- 
tion and heating associated with mass transport occurs 
in an extremely non-local fashion. While it seems im- 
possible to avoid some critical /3, its actual value need 
not be near one, but rather is given by Equation ([TQ|) . 
modulo small modifications due to spiral structure. 

It is dangerous to draw conclusions from un-converged 
simulations, but if the conversion between Gl-driven tur- 
bulence and heating is treated properly in these runs, we 
suggest that rather than calling the Viscous Criterion 
into question, it implies that a^at is decreasing with in- 
creasing resolution, (a possibility that remains untested). 
One should not be surprised if Ofgat remains un-converged 
at the maximum resolution, as the effective resolution is 
only marginally higher {5x = 1/ 3H) than that at which 
c onvergence was ac h ieved in 2D ()Gammiell2001f ). 

iLodato fc Clarkd (|2011l ) have suggested that devia- 
tions in a due explicitly to artificial viscosity are only a 
few percent. However, it seems plausible that the resolu- 
tion is simply too low to resolve the turbulence properly 
in the first place (another possibility which they acknowl- 
edge). If the highest resolution calculations are converg- 
ing on an accurate value of agat , it must be understood 
why it is a fact or of 10 — 20 smaller than that found by 
iGammid ()2001f ). The dimensionality of the calculations 
could explain the discrepancy; similar critical /3 values 
were found in 2D by iGammid (^2001) and in low resolu- 
tion SPH simulations which a re effectively 2D due to a 
lack of scale height resolution (jRice et al.ll2003f ). 

If the high resolution simulations have low o^sat's that 
are representative of the maximum transport provided by 
GI in low mass, thin disks, then this is worthy of further 
investigation. For sufficiently small agat, the Stalling Cri- 
terion may take over as the dominant cooling criterion. 

8. SUMMARY AND DISCUSSION 

In this paper, we separate the conditions for successful 
fragmentation into two considerations: (1) Can the disk 
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reach the threshold of hnear instabihty? and (2) Can a 
fragment survive after it has begun to cohapse? This is 
distinct from previous work that has focused on whether 
or not the disk can reach Q = 1 and cool quickly. We 
show that cooling is intimately connected both to reach- 
ing instability and to fragment survival. 

We first address the distinction between disks heated 
internally by the dissipation of accretion energy, and 
those heated externally by irradiation. We show that 
the Viscous Criterion only governs the answer to ques- 
tion (1) in viscously heated disks because irradiated disks 
always effectively satisfy the Viscous Criterion. We show 
that fragmentation is easier in irradiated disks, and that 
if an irradiated disk is thin enough (or equivalently low 
enough in mass) when Q ~ 1, global transport may be 
ineffective, allowing the disk to be driven to collapse by 
low rates of mass infall. Thus, fragmentation may occur 
at slightly smaller radii, and lower initial masses then 
previously assumed. Importantly, fragmentation at low 
accretion rates could reduce the expected growth of frag- 
ments after their formation. 

We then derive a generic cooling criterion which con- 
trols fragment survival in all disks. The Stalling Crite- 
rion divides fragments that can collapse on their own 
free-fall timescale from those that will become pres- 
sure supported and begin Kelvin-Helmholtz contracting. 
Fragments in the latter category must contract to a criti- 
cal radius (i?H+p) to avoid being torn apart by tides when 
pressure stalls their collapse. This radius is roughly half 
the Hill radius: the reduction is due to the inclusion of 
pressure support in the calculation of force balance. 

For reasonable perturbations, fragments in disks with 
soft equations of state (7 = 7/5) will likely be able to 
contract below i^H+p, even if they cool over many or- 
bital periods. In contrast, for stiffer equations of state, 
fragments that do not collapse in free fall will likely be 
sheared apart. 

We note that the value of Ru-\-p depends on fragment 
geometry, and whether a given fragment will contract 
to this size depends on its non-linear evolution in time. 
For example, if fragmentation occurs at densities and 
temperatures where there are sharp jumps in opacity, 
this must be taken into account in deter ming /3. Specific 
contraction models for disk-born fragments at all masses 
should be compared to the Stalling Criterion. We com- 
ment on other disruptive processes, including the classic 
Collisional Criterion, and the effects of fragment rota- 
tion, in Section [5l 

Of course, since fragments will continue to grow fol- 
lowing formation, gravitational instability should still fa- 
vor the production of low-mass brown dwarfs over plan- 
ets (Kratter et al. 2010b); That GI produces larger ob- 
jects is increasingly supported by observati ons of M-star 
and brown- dwarf companions to A-stars ( Hinklev et al.l 
I2OIOD, and the paucity of wide, planetary mass compan- 
ions (jLeconte et al.ll2010f ). 

8.1. Suggestions for Future Work 



This work should be a useful guide for future nu- 
merical simulations intended to study fragmentation in 
irradiation-dominated disks. Two major uncertainties 
limit our conclusions for protostellar disks: 

1. What is the value of agat in a viscously- heated Q ^ 
1 disk? 

2. What is the maximum Gl-driven transport rate as 
a function of H/r in an irradiated Q ^ 1 disk? 

These two uncertainties are effectively the same question, 
but since the physics of transport in these regions differs, 
they must be considered independently. Simulations will 
also be required to confirm our conclusion that regard- 
less of whether they collapse in free fall, fragments can 
likely survive if 7 = 7/5 because they collapse to smaller 
than i?H+p on of order a free-fall time. Hence, disks 
with soft equations of state must simply reach Q = Qo 
to fragment. Confirmation of this conclusion will require 
a full treatment of radiative transfer and disk opacity. 
As noted earlier, the impact of fragment rotation must 
also be explored. Current state of the art simulations of 
disks with radiative transfer do not yet cover the range of 
parameter space required to answer the above questions 
(Boss 2007; Cai et al. 2008, 2010; Boley et al. 2010; note 
that the latter two authors have questioned the conclu- 
sions of the former based on differing treatments of the 
photosphere). 

Answers from such simulations will enable us to dis- 
tinguish between the following outcomes in each regime. 

• Viscously-heated regions: If a^at ^ O-l, then stan- 
dard applications of the cooling time criterion give 
accurate results. If asat ^ O-l^ then the standard 
Viscous Criterion will determine the boundary for 
fragmentation, but the critical cooling time will be 
longer. 

• Irradiation- dominated regions: If both the maxi- 
mum Gl-driven transport rate and aeff due to MRI 
are too low to process the infall rate onto a pro- 
tostellar disk, fragmentation may be possible for 
thin disks at infall rates lower than ^cl/G. If this 
occurs, low-mass fragments which do not experi- 
ence substantial growth might be produced. Oth- 
erwise, high accretion rates are required to drive 
these disks unstable. 
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APPENDIX 
BEHAVIOR OF Q WITH M 

To determine under what conditions increasing M drives disks more unstable, we consider the time evolution of M, 
Q, and T. Reformulating Equation (|4]) , the disk accretion rate becomes: 

M = 3^§, (Al) 

where we have assumed that a is some unknown function of Q f Lin fc PringldflQSTl : iKratter et alJ[2QQ8[ ). Because the 
evolution of Q depends on how accretional heating is converted into midplane temperature, which is governed by the 
disk opacity, we must consider a temperature dependent opacity: 

K = noT^-. (A2) 

Solving for Q/a{Q) and using Equation (J5j) to get T as a function of M we have: 

10-2/3^ 12-3/3,^ 

Q f3\ '-'^- fkB\ '-^^- f 1^0^^ , - I - \ - I - \ /AQ^ 

[A6) 




a \G J V M / \167r^(j 



where we have rewritten S using the definition of Q. We subsume the first three terms of Equation ()A3p into a constant 
ai for convenience. Rearranging, we find: 

Q = af^a^^^M^ (A4) 

Because a is an unknown function of Q, we relate the time derivatives of a and Q using an unknown function Ci{Q) > 0: 



-^ = -cA^. (A5) 

a dt Q dt ^ ^ 

This relation captures the fact that we expect a to increase as Q declines, and allows us to treat the above expression 
analytically. Taking the time derivative of Equation ()A4p we find: 

l0-2p^ ^) Q dt b-p^M dt' ^ ^ 

We can immediately see that the sign of the relation between the derivative of M and the derivative of Q is governed 
by the value of P^^ and in some cases Ci. To aid in our interpretation of Equation ()A6p . we can also relate dT/dt and 
dM/dt: 

T dt \ 5-f3^Ci{7/2-f3)J 5-f3M dt ^ ^ 

Together these two equations dehneate four different f3^ regimes: 

1. /3^ < 2: In this case, the right hand side of Equation (|A6[) is negative, so that higher accretion rates drive the 
disk more unstable. Thus for small /3,^, an increase in M causes a modest increase in temperature and a larger- 
increase in S, causing Q to decline. This case is most relevant to protostellar disks. See iMatzner fc LevinI (|2QQ5l ) 
for a discussion of the P^ = 2 case. 

2. 2 < /3,^ < 7/2: Both sides of Equation ()A6p are positive, meaning that an increase in M actually drives Q up. 
This occurs because the steep increase in opacity with T causes the temperature to rise faster than H, increasing 
Q. 

3. P^ > 5: When /3,^ is extremely high, an increase in M always corresponds to a decrease in Q. This dependence 
derives from the fact that for /3,^ > 5, an equilibrium between accretional heating and cooling as defined by 
Equation ([5]) requires that an increase in M correspond to a decrease in temperature, which is counterintuitive 
(T (X Ms^^ at constant a). 

4. 7/2 < Pf^ < 5: The sign of the left hand side of Equation (|A6[) depends on the value of Ci , and thus the functional 
form of a{Q). Incr easin g Ci shifts the value of /3 at which an increase in M corresponds to a decrease in T, as 
shown in Equation (|A7[) . When Ci ^ 0, the disk behaves as in case 2 all the way to /3 = 5. When Ci is large, 
this case behaves as case 3, with the upper limit on /3 a function of Ci. 

The outer regions of protostellar disks on which we focus here are all encompassed by case 1, in which increased 
accretion destabilizes the disk. 
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